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Space-time dependence of the anomalous exponent of electric transport in the disorder model is 
presented. We show that the anomalous exponent depends on time, according to the time-evolution 
of the number of the effective neighbouring sites. Transition from subdiffusive to normal transport 
at long-enough time is recovered. The above result indicates that the spatial structure, specifically 
the network structure, of the hopping sites might be a novel element which determines the anomalous 
exponent. This leads to the feature that the scaling property of the electric transport to time is 
insensitive to other elements, such as the distance of the sites or the spatial dimension of the system. 
These findings are verified by means of Monte Carlo simulation. The relation of the result to the 
conventional knowledge of the Multiple Trapping Model is shown by deriving it as a special case of 
the disorder model. 

PACS numbers: 72.80.Ng, 02.50.-r, 05.40.Fb 
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INTRODUCTION 

Noncrystalline materials, representative examples of 
which are organic electro-luminescence, organic pho- 
toconductors, and amorphous Si photoconductors, are 
readily seen around us. For instance, photoconductors 
are used in photocopiers, laser beam printers, and solar 
cells. Since to control the electric conductivity in the 
above-mentioned devices is crucial for their development 
and design, it is necessary to investigate the mechanisms 
of their carrier transport. 

Electric conductivity in noncrystalline materials is 
known to be dominated by the so-called hopping conduc- 
tance. It is widely recognized that hopping conductance 
results in anomalous diffusion (subdifFusion) at a certain 
mesoscopic scale [1]. Experimentally, the anomaly man- 
ifests itself in the long-tail of the time-of-flight (TOF) 
signal [1-3]. 

Anomalous diffusion has been thoroughly studied in 
the context of continuous-time random walk (CTRW) [4] . 
In CTRW, the waiting time of the walker depends 
on time as 'w{t) ~ in the limit t oo, where 

a is called the " anomalous exponent" . SubdifFusion cor- 
responds to the case < a < 1, while a > 1 describes 
superdifFusion, or the Levy flight. It is known that the 
continuum limit of CTRW is described by the fractional 
diffusion advection equation of order a [4] . The crucial 
feature of CTRW is that the above-mentioned anomaly 
can be expressed by a single parameter a. In fact, by 
varying a, the long-tail of the TOF signals can be fitted, 
as shown in ref. [1] . However, since a is treated as a mere 
parameter in the CTRW, it is significant to study the re- 
lation of a and physical quantities of interest, such as the 
density of state (DOS) of the trap levels, or the spatial 
structure of the hopping sites, for physical understanding 
of the hopping conductance. 



Almost 30 years ago, such a relation was partially given 
by the "multiple trapping model (MTM)" . In this model, 
a is related to the DOS of the trap levels and the tem- 
perature T as a = T/Tc, where Tc is the typical width 
of the DOS. The DOS is assumed to be of the exponen- 
tial type, p{E) = e^/C^^^^VlfeB^c) {E < 0), where the 
origin of the energy is set to the edge of the conduction 
band. This type of DOS is typical for disordered irror- 
ganic semiconductors [5]. 

On the other hand, hopping conductance has been 
studied mainly by two models, the "disorder model [6]", 
which is studied in this paper, and the "polaron model 
[7]". The former is suitable for a system where the 
anomalous charge transport is dominantly activated by 
the static energy disorder of the hopping sites, and where 
the effect of the weak electron-phonon coupling is rela- 
tively negligible. In contrast, the latter is suitable for a 
system with strong electron-phonon couplings, and with 
relatively negligible effects of the energy disorder [8] . 

In the disorder model, the hopping rate between two 
sites is described by the classical model of [9], which is 
based on [10]. As we will review briefiy, this model in- 
corporates the information of the spatial structure of the 
hopping sites, which the MTM does not take into ac- 
count. To solve the model of [9], Monte Carlo (MC) 
simulation has been performed [6], and the relation be- 
tween the disorder model and the MTM has been studied 
numerically [11]. However, the relation of the anoma- 
lous exponent and the physical quantities, including the 
spatial structure of the hopping sites, has not been es- 
tablished theoretically until present. Noting that the re- 
lation of the anomalous exponent and the width of the 
DOS has already been established in the MTM, the the- 
oretical understanding of the "disorder model" seems to 
be left in a premature state. The aim of this paper is to 
promote this theoretical understanding. 

In this paper, we first review the disorder model and 



2 



its related issues. Secondly, we show the premise and 
the main result with emphasis on the novel aspects of 
our study. Then, the theoretical result is verified by MC 
simulation. The methodology of the simulation is briefly 
explained and the consistency of theory and simulation 
is exhibited. Next, the relation of our result and that of 
the MTM is revealed. Finally, we summarize our study. 



THEORY 

In this study, we consider the disorder model. We 
adopt as the DOS the exponential type mentioned above, 
in order for the correspondence to the MTM. For other 
types of DOS, e.g. Gaussian DOS, which is typical for 
organic semiconductors, similar techniques would be ap- 
plicable. 

Before we proceed to the derivation of the anomalous 
exponent, we briefly review the waiting time w{t) of the 
CTRW and the hopping probability of the model of [9]. 
In CTRW, the probability distribution that the carrier is 
at position x at time t is given as 

p{x,t) = y2 / dTp{x',T)il){x,x',t-T) + ^{x,t)5xo, 

l-L Jo 

if the carrier is situated at a; = at f = 0. Here, ^{x, t) = 
1 — Jq dTw{x,T), and %l}{x',x,t) denotes the probability 
density that the carrier hops from x to x' after waiting 
time t. The waiting-time probability distribution w{x,t) 
is given in terms of ijj{x' ,x,t) as 



wix,t) = ^ i){x',x,t). 



For the disorder model, ip{x',x,t) is expressed in the 
splitted form, tp{x',x,t) = 'w{t)(p{x' ,x), where 4){x' ,x) 
is the spatial part of the probability distribution, which 
is normalized as X^^,/ x) = 1. The time-dependent 
part w{t) is of our interest. 

The hopping probability, i.e., the probability per unit 
time of the carrier to hop from site i to site j, which is 
denoted Vij, is given as [9] 



where = 'Ylij^j^'^ij- Here, N represents the set of 
all the sites in the system. Eq. (2) indicates that the 
hopping process is modeled as a homogeneous Poisson 
process with decay rate A^. The physical origin of the 
disorder of energy is that of the spatial distribution of the 
hopping sites. In the disorder model, the spatial disorder 
is assumed to be moderate enough that the positions of 
the sites fluctuate around the structured lattice points, 
according to some distribution (e.g., Gaussian distribu- 
tion). However, for the sake of simplicity, we incorporate 
the effect of the spatial disorder to that of the disorder 
of the site energy, and assume that the sites compose a 
structured lattice with lattice spacing a. 

Now we calculate the anomalous exponent. We con- 
sider a system with arbitrary spatial dimensions, where 
the carrier at each site is allowed to hop to all the sites in 
the system. The crucial observation is that the waiting- 
time distribution Eq. (2) is merely valid for a certain en- 
ergy. If the energy is distributed according to some prob- 
ability distribution, then the ensemble average with re- 
spect to the energy distribution, which we denote {w{t)), 
should be regarded as the effective waiting time of the 
system [12]. 



{w{t)) = n / de^jPLie^j) Wi{t) 
= n / de,' 



2kBTc 



X Y^KkiB e-^'^^^"^ * , (3) 



where 



(4) 



, and eji = Ej — Ei is the difference of the energy of site i 
and j. Note that translational invariance is assumed, as 
a result of the spatial coarse-graining due to the integra- 
tion with respect to site energies, and the siibscript i is 
omitted in the expression {wit)). In Eq. (3), we have uti- 



v^j{R,j,Ej - Ei) = i,„e~'^^^*i-('^i-^i)®(^i-Ei)/kBTi^-^^^ lj2;ed the fact that the energy difference ej^, which is the 

difference of two identical probabilistic variables, obeys 
the Laplace distribution PL{eji) = e'^'^^^/''^^'' /{2kBTc). 
We can calculate the integrals in Eq. (3) by two suc- 
cessive transformations of the integration variables, i.e., 
Aji = e~'^i'/^'^'^ for the first step, and Cji = KjiAjit 
for the second [13]. By the first transformation Aji = 
g-eji/ksT^ and using dAji/dcji = —Aji/ksT, we obtain 



Here, Ei is the energy of site i, Rij is the distance of site 
i and j, £, is the damping factor of the wave function in 
the localized state, vq is a coefficient which represents the 
magnitude of the hopping rate (typically of the order of 
10^^[sec~^]), 9(a;) is the Heaviside's step function, and T 
is the temperature. Dependence on microscopic physics, 
such as the carrier-phonon interaction and the phonon 
DOS, are assumed to be incorporated in i/q- The waiting 
time of the carrier on site i, which we denote Wi{t), is 
expressed in terms of Vij as 

Wi{t) = Aie-^*\ (2) 
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Then, by the second transformation Cji = KjiAjit, we 
obtain 



1 T 



Kut 



IT T ( T 

= 2 7r(^iii)"'^7(^,^.,i 



Here, 



^{T/T,,Kijt)= / drr-i 
Jo 



(6) 



(7) 



is the lower incomplete gamma function. Similarly, we 
obtain 



Jo 



> ,, ^ p-\<iji\/kBTc 



(2fcBTe) 



Substituting Eqs. (6) and (8) into Eq. (3), we have 



xn 



jeMdi^k 2 



T 7(^ + l,-ft^ifc*) „-Kikt 
Tc {Kikty+T/Tc "I" ^ 



Tc (Kijt)T/Ta 



(8) 



(9) 



It is important that the following relations hold for the 
sites which satisfy the condition Kikt <C 1, 



T 7 



kt 



1, 



+ e 



1. 
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Furthermore, for the sites which also satisfy the condition 
Kik -C 1, we obtain 



2 



T 7 + l^-fi^jfci 



0, (12) 



in combination with Eq. (11). Therefore, the sites which 
satisfy ii'ifet ^ 1 and Kik ^ 1 arc irrelevant to the effec- 
tive waiting time. Then we can rewrite Eq. (9) as 



K, 



ik 



n I 



y 7(^^ + 1, K.fei 

{K,kty+'^'/^^ 



+ e" 



-Kikt 



+ e- 



(13) 



Here, AfR{t) is the set of the "effective neighbours" at 

time t, which is defined by removing the sites which sat- 
isfy K^kt < 1 and K^k < 1 from Af. 



Next, to obtain further an approximate analytic ex- 
pression for the effective waiting time, we attempt to ap- 
proximate the lower incomplete gamma function by the 
gamma function. Prom Eq. (7), it is obvious that 



j{T/n,Kijt)c^T{T/T,) 



(14) 



for Kijt 3> 1. To be more specific, we introduce a crite- 
ria which defines the precision of the approximation Eq. 
(14). We consider a constant Cj, where the following 
replacement is performed for the CcLS6 Kijt C/yl 

7 (T/Te, Kijt) ^ r (T/Te) {Kijt > Cj,0 < T/T, < 2). 

In other words, defines the precision of the approxi- 
mation Eq. (14). Then, we can define the critical radius 
R{t) at time t, 



R{t) = - 



ln% 
2^ ■ 



(15) 



For sites with Rij < R{t), approximation Eq. (14) is 
applied in Eq. (13). 

We approximate Eq. (13) by only including the neigh- 
bouring sites with Kijt > Cj. Then, 



'kit) 



2 



T \ r( j^) -Kikt 

1 + T/Tc 



Taj (K,kt) 

T r(j,) 



-Kiit 



(16) 



where A/'^(t) is defined as a set of neighbouring sites with 
Rij < R{t). Let us denote the set of the n th- nearest 
neighbouring sites as A/"™. From the definition of Kij , Eq. 
(4), it is obvious that Kikt = e'^^"^°-Kijt (k e Afn-m, j S 
N'n, m eN) holds for the n th-nearest and the {n—m) th- 
nearest neighbours. For a system with typical parameters 
(e.g., the case shown in the numerical simulation later), 
g2{a ^ ^ holds. Hence, for time t where the n th-nearest 
neighbours are excluded in Eq. (13), i.e., Kijt ^ 1 (j e 
Mn), there exists m G N such that Kijt > Cj (j S 
J\fn-m), i-e., {n—m) th-nearest neighbours are included in 
A/]j(t). For the case m= 1, {n — 1) th-nearest neighbours 
arc included in Af^{t), and hence A/ij(f) in Eq. (13) and 
M'fi{t) in Eq. (16) coincide. In general, the (n — m) th- 
nearest neighbours with 1 < Kijt < Cj is omitted in Eq. 
(16), but as far as e^^" 3> 1 holds, their effect on the 
effective waiting time is subdominant compared to that 
of the sites with Rij < R{t). 

We can introduce a characteristic time t„, at which 
the n th-nearest neighbours contribute to the efi'ective 
waiting time, as 



Tn = 



-2^an ' 



(17) 



From Eq. (17), we obtain t„+i = e^"T„. As mentioned 

above, r„+i ^ r„ holds for typical systems. This indi- 
cates that the {n+1) th-nearest neighbours are negligible 
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for time scales around or below r„. Thus, we can define 
a time interval Jjv where the number of the "effective 
neighbors" equals N. Then, it can be stated that Eq. 
(16) is valid for t <E JjVj,(t), where Nfi{t) is the number 
of sites in M^{t). From Eq. (15), the critical radius R{t) 
increases monotonically with time. As a result, Nfi{t) 
increases monotonically as limt^oo N'ji{t) = Ntot, where 
Ntot is the total number of sites in the system. 

Prom Eq. (16), we can derive simple results for the fol- 
lowing two limiting cases: NR{t)T < and Nn,{t)T > 
Tc. For the case N[{{t)T > T^, taking the limit T ^ oo 
with Tc fixed in Eq. (16) leads to 



{w{t)) 



1 



' 2^36Ar' (t) ^i]*- 



(18) 



This is the usual Poisson distribution, and hence a = 1, 
i.e., the diffusion is normal. This is consistent with our 
understanding that the band conduction dominates at 
high temperatures. 

For the case Nii(t)T <C Tc, one can see that the first 
term is dominant, respectively, in the two squared brack- 
ets in Eq. (16): 

xn,w;.(*).^.(^^r(^)(i^,,i)-^)} (19) 



oc 



The asymptotic form (wit)) oc i-(i+WR(*)T/T,) ^ 
7jVH(t)) indicates that the coarse-grained hopping conduc- 
tance possesses the feature of subdiffusion with a time- 
dependent anomalous exponent a(t) ~ Nji{t)T/Tc. 

To summarize, we derived the asymptotic form of the 
anomalous exponent. 



a{t) = 



1 ° (Arfl(t)T»Te) 



(20) 



which is valid for time t G lNR{t)- Eq. (20) indicates 

that the anomalous exponent does depend on the spa- 
tial structure of the hopping sites. Moreover, an inter- 
esting feature of Eq. (20) is that the network structure 
of the hopping sites is the only spatial element which 
contributes to the anomalous exponent. In other words, 
other elements, such as the distance between the sites 
and the spatial dimension of the system, are insensitive 
to the time-dependence of the electric transport, as long 
as Nii{t) is identical. It is also clear from Eq. (20) that 
the anomalous exponent depends on time. Since Nu{t) 
increases with time, a{t) eventually becomes 1 at some 
time in the future, i.e., subdiffusive transport necessar- 
ily becomes normal. This is consistent with the common 
knowledge of nonequilibrium statistical mechanics that 
the Markovian property is recovered by a proper coarse 
graining [14]. A novel indication of the above result is 
that, since Nfi{t) increases with time, a discrepancy with 



the MTM, which corresponds to the case of Nji{t) = 1, 
shows up. This issue will be discussed later. 



SIMULATION 

The theoretical asymptotic form of the anomalous ex- 
ponent Eq. (20) is verified by the MC simulation of the 
hopping conductance. We mainly focus to the following 
three theoretical predictions : (i) the transition from the 
subdiffusive transport to the normal transport, (ii) the 
transient behavior of the anomalous exponent, and (iii) 
the effect of the network structure of the hopping sites 
to the anomalous exponent. 

The anomalous exponent results in the time- 
dependence of the transport coefficients. This fact 
suggests that one cannot adopt the naive defi- 
nition of the transport coefficients, e.g., D = 

[{x'^(ttot)) — {x{ttot))'^] ^ /ttot for the diffusion coefficient 
[6], where ttot is the total duration time of the hopping 
process. Rather, one should estimate the corresponding 
moments at intermediate times, e.g., 

D{t) = [{x\t)) - {x{t)f] yt {0<'t< ttot). (21) 

Estimate such as Eq. (21) cannot be performed in the 
conventional MC methods [6], which is difficult to syn- 
chronize the hopping procedure of the independent car- 
riers. An algorithm is developed to synchronize the hop- 
ping of the carriers, whose details will be published else- 
where [15]. We adopted this algorithm to calculate the 
time-dependence of the transport coefficients, or, cquiv- 
alently, the anomalous exponent. To be more precise, in 
order to estimate a from the MC simulation results, we 
used the following relation which holds for subdiffusion: 



{x^it))-{x(t)f io<yt<ttot). 



(22) 



From Eq. (22), we can estimate a by fitting the function 
f{t) = Ft", where F is a constant, to the simulation re- 
sult at time t. This estimate requires the synchronization 
of the hopping carriers. 

Before proceeding to discuss the simulation results, we 
estimate the value of in Eq. (15). Since the incom- 
plete gamma function j {T/T(.,C^) asymptotically ap- 
proaches to the gamma function P (T/T^) , at least in 



6-7 digits accuracy, at = 100, we chose C-y 



100 



for the case 0.05 < T/Tc < 1. Substituting typical 
values t = 10-2 [sec], 2^a = 10, a = 10"^ [m] and 
z/o = 10"'^^ [sec"-'^] to Eq. (15), we obtain the critical 
radius i? = 2.8 x lO'^ [m] for = 100. This indicates 
that it is enough to consider at most the third-nearest 
neighbours for typical systems. Therefore, in realistic 
time scales, the number of the effective neighbours is fi- 
nite, typically at most of the order of ten. This is consis- 
tent with common hopping models, which assume that 
the number of the neighbouring sites to which hopping is 
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allowed is finite [6]. In addition, using the above param- 
eter values, we obtain a relation of characteristic times, 
Tn+i — e'^^°-Tn — 2.2 X 10^T„. As mentioned in the pre- 
vious section, this confirms that {n+ l)th-nearest neigh- 
bours are negligible for time scales around or below r„. 

Now, we verify the validation of Eq. (20) by the 
MC simulation. First, in order to illuminate the tran- 
sition from subdifFusive to normal transport, we chose 
a one-dimensional system with carriers allowed to hop 
up to the second-nearest neighbouring sites. The condi- 
tion of the simulation is as follows; the number of car- 
riers is Np — 1000000, the parameters are 2^a = 10, 
I/O = 10^^[sec-i], T/Tc = 0.375. These values corre- 
spond to the characteristic times ri = 2.2 x 10~^[sec], 
T2 — 4.9 X 10^^ [sec], etc. Initially, all the carriers are 
rested at the origin. For the case where the correlation 
between the carriers can be neglected, this initial con- 
figuration is sensible, since in the novel algorithm [15] 
neighbours are generated independently for each carrier, 
where each neighbour spans a probabilistically indepen- 
dent sample space. 

We can theoretically estimate characteristic values of 
a(t) at characteristic times. Since T/Tc = 0.375 and 
there are two sites in the 1st and 2nd-nearest neighbours 
in the system, respectively, a(t) becomes as follows: 



ait) 



0.75 (t = Ti) 

1.0 {t>T2) 



(23) 



Eqs. (23) indicate that a{t) approaches 1.0 eventually, 
taking the value 0.75 at t = ri. 

The result of the MC simulation is shown in Fig. 1. 
One can see that the result is compatible with the theo- 
retical estimate, Eqs. (23). Thus, transition from subd- 
iffusive to normal transport is verified. Moreover, it can 
be seen that the transport is normal in the early stage of 
the diffusion. In this stage, since only carriers with short 
life-times can hop to the neighbouring sites, the hopping 
is approximately unimodal. Therefore, the transport is 
almost normal. 
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FIG. 1. The result for the MC simulation of the transition 
from subdiffusive to normal transport. The system is one- 
dimensional with hopping up to the second-nearest neigh- 
bours is allowed. 

Next, in order to illuminate the transient behavior of 
the anomalous exponent, we chose a one-dimensional sys- 



tem with carriers allowed to hop up to the third-nearest 
neighbouring sites. The condition of the simulation is 
as follows; the number of carriers is Np ~ 1000000, the 
parameters are 2^a = 10, = lO^^jggc-i], T/T^ = 0.1. 
These values correspond to the characteristic times ri = 
2.2 X 10"S[scc], T2 = 4.9 X 10"2[sec], = 1.1 x 103[sec], 
etc. 

Characteristic values of a{t) at characteristic times can 
be estimated theoretically. Since T/Tc — 0.1 and there 
are two sites in the 1st, 2nd and 3rd-nearest neighbours 
in the system, respectively, a(t) becomes as follows: 



ait) 



0.2 (t^Ti) 

0.4 {t = T2) 

0.6 {t = T3) 



(24) 



The result of the MC simulation is shown in Fig. 2. 
One can see that the result is compatible with the the- 
oretical estimate, Eqs. (24). In addition, plateaus are 
found around each characteristic time r„ {n = 1,2,3). 
This indicates that a{t) can be regarded as constant at a 
certain time interval. This is compatible with the com- 
mon treatment that a{t) is regarded as constant [1]. 
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FIG. 2. The result for the MC simulation of the transient 
behavior of the anomalous exponent. The system is one- 
dimensional with hopping up to the third-nearest neighbours 
is allowed. 

Finally, in order to illuminate the effect of the network 
structure of the hopping sites, two cases are chosen for 
illustration. One is a one-dimensional system with car- 
riers allowed to hop up to the second-nearest neighbours 
(4 sites in total, see the upper figure in Fig. 3). Another 
is a two-dimensional system with carriers allowed to hop 
only to the nearest neighbours (4 sites in total, see the 
lower figure in Fig. 3). The condition of the simulation 
is as follows; the number of carriers is Np = 1000000, 
the parameters are = I, i^o — 10^^[sec~^], and T/Tc 
is varied in the range 0.025 - 1.0. The duration of the 
simulation is chosen to be larger than T2. From Eq. (20), 
it is expected theoretically that the results for the two 
cases are identical. 

The result of the MC simulation is shown in Fig. 3. 
One can see that the results for the two cases are identi- 
cal. For both cases, the simulation result coincides with 
the theoretically derived line a = AT/Tc for T/Tc < 0.2, 
and it reaches a = 1 for T/Tc > 0.6. The network struc- 
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tures of the two systems are identical, i.e., they are home- 
omorphic, while other features such as the distance be- 
tween the sites and the spatial dimension differ. This 
result indicates that the anomalous exponent is deter- 
mined by the network structure of the hopping sites. 
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FIG. 3. Comparison of the theoretical and the simulational 
results for the relation between the anomalous exponent a 
and T/Tc- The upper figure is for a one-dimensional system 
with hopping up to the second-nearest neighbours are allowed 
{Nu{oo) = 4), while the lower figure is for a two-dimensional 
system with that only to the nearest neighbours are allowed 
{Nu{oo) — 4). The network structure of the two systems are 
identical (homeomorphic) , while other features such as the 
distance of the neighbouring sites and the spatial dimension 
of the system are different. 

The origin of this network structure dependence re- 
sides in the feature of the model that the probability 
distribution of the waiting time is described by a Pois- 
son distribution Eq. (2) for a given site energy. Due to 
this feature, the dependence of the time-dependent part 
of {'w{t)) on the distance of the sites appear (through 
Kijt) in the upper limit of the integral of the incom- 
plete gamma function, which is wiped out to infinity in 
the long-time limit, leaving the power-law dependence of 
Kijt. The emergence of the Poisson distribution is at- 
tributed to the independence of each hopping procedure, 
which is valid in the weak-coupling limit of the carri- 
ers. For the strong-coupling case, intricate distribution 
will possibly emerge, and it is beyond the scope of the 
present paper. 



DISCUSSION 

One of the major results of this paper is that the 
anomalous exponent in the disorder model [9] is in gen- 
eral larger than that in the MTM [16]. To the extent of 
our knowledge, an exact correspondence between the dis- 



order model and the MTM has not been stated clearly. 
In the remainder we will state their correspondence. 

In the MTM, the carriers reside in the infinite tower 
of almost-continuous energy levels ("trapped levels"), as 
well as in the conduction band. The carriers in the con- 
duction band ("conduction carriers") and those in the 
trapped level ("trapped carriers") are treated separately. 
The "conduction carriers" are allowed to move to the 
neighbouring sites according to normal diffusion, while 
there is no diffusion for the "trapped carriers" . 

On the other hand, there is no distinction between the 
"conduction carriers" and the "trapped carriers" in the 
disorder model. In particular, the "trapped carriers" are 
allowed to hop to the neighbouring sites as well. 

These two models apparently seem independent, but 
by imposing a constraint on the disorder model that hop- 
ping to the neighbouring sites are allowed only when the 
energy level of the sites coincide with (the energy of 
the lower edge of the conduction band), the anomalous 
exponent of the MTM can be obtained. To illuminate 
this correspondence, we chose a one-dimensional system 
with a carrier allowed to hop only to the nearest neigh- 
bouring sites. Imposing the above constraint to Eq. (1), 
the hopping rate from site i to site i±l becomes 

,-(E^^E,)/kBT 





Vie 



Ec) 



{Ei±i < Ec) ^ ^ 



Here, Vi = VQe~^^°'. Using Eq. (25) and choosing Ec = 
for simplicity, we obtain the waiting time for this case w' 
as follows: 



w'{E,^uE„E,+i) 

2y^^E^lksT^-2.^e^^l^BT, ^^^^^ ^ ^ 



(.26) 







(S,+i < 0, E,^i < 0) 



By averaging over the energy, we can obtain the effective 
waiting time as 



{w'{E,^i,E„E,+i)) 







i-l 







ksTc^ 







ksTc 
Ei+i 
ksTc 



(27) 
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Here, 0_ = — lima^+o ^- Note that, due to the intro- 
duction of a specific energy Ec{= 0), the integral in Eq. 
(3) should be performed independently for the three sites 

i,i ± 1, in the region Ei,Ei±i G (— oo,i?c), respectively. 
Now, taking the Imit t ti, and using Eq. (8), we have 

+2.,^M)---rg + l)lin.e^ 
(xt-'^-T^. (28) 

From Eq. (28), wc have a = T/Tc- This is exactly 
the result of the MTM [16, 17]. The generalization of 
the above result to the system with arbitrary numbers 
of neighbouring sites is tedious but straightforward , and 
the same result can be obtained. Thus, the MTM can be 
expressed as a special case of the disorder model. Prom 
the viewpoint of the disorder model, the spatial factor 
Nji{t) in the anomalous exponent is exactly cancelled by 
the constraints imposed above, leaving a to be superfi- 
cially insensitive to the spatial structure. 



SUMMARY 

In this study, we have calculated the anomalous ex- 
ponent in the disorder model. First, we have shown 
that the anomalous exponent is time-dependent. Fur- 
thermore, since the characteristic time scales t„ where 
the n th-nearest neighbours contribute to the effective 
waiting time are well separated, i.e., r„+i ^> r„, the 
{n + 1) th-nearest neighbours can be ignored for time 
scales around or below r„. This indicates that the num- 
ber of the effective neighbouring sites Nfi (t) is finite in re- 
alistic (experimental) time scales. In fact, we have shown 
that it is enough to consider at most the third-nearest 
neighbours for typical systems. On the other hand, since 
Nji{t) increases monotonically with time, a{t) eventually 
becomes 1 at some time in the future, i.e., subdiffusive 
transport eventually becomes normal. This is consistent 
with the common knowledge of nonequilibrium statistical 
mechanics. 

Secondly, we have shown that the spatial structure of 
the sites might be a novel and an equally significant ele- 
ment in the disorder model, in addition to the two com- 
mon elements, the temperature T and the typical width 
of the DOS, Tg. Specifically, the network structure of 
the sites is shown to be the novel element, which can 
be seen in the expression for the anomalous exponent, 
a{t) = Nii{t)T/Tc, where NR{t) is the number of the 
"effective neighbouring sites" at time t. This leads to 
the feature that other elements, such as the distance be- 
tween the sites or the spatial dimension of the system, are 



insensitive to the time-dependence of the electric trans- 
port, as long as iVij(f) is identical. 

We have verified the above theoretical results by means 
of MC simulation of the hopping conductance. First, we 
chose a one-dimensional system with carriers allowed to 
hop up to the second-nearest neighbouring sites, to verify 
the transition form subdiffusive to normal transport. The 
result actually showed the transition, and the anomalous 
exponent was identical to the theoretical values. Next, 
we chose a one-dimensional system with carriers allowed 
to hop up to the third- nearest neighbouring sites, to ver- 
ify that (n + l) th-nearest neighbours can be ignored for 
time scales around or below t„ . The result was consistent 
with the theoretical prediction. Moreover, plateaus were 
found around each r„. This indicates that the anomalous 
exponent can be regard as constant for certain time in- 
tervals. This is compatible with the common treatment 
where the anomalous exponent is assumed to be constant. 
Finally, a one-dimensional and a two-dimensional system, 
with identical network structures, were compared, to ver- 
ify that the network structure determines the anomalous 
exponent. The results were identical for the two cases, in 
agreement with the asymptotic results for the anomalous 
exponent, Eq. (20). 

We have also revealed the relation between the disorder 
model and the MTM [16, 17], which was not stated and 
remained obscure before. We have shown that the MTM 
can be expressed as a constrained version of the disor- 
der model, where only hopping to the energy levels of 
the conduction band is allowed. This constraint exactly 
cancels the factor iV/j(i) in the expression of the anoma- 
lous exponent. Whether the system is better described 
by the "general (unconstrained) " disorder model or the 
"constrained" disorder model (i.e., MTM) is a matter of 
the nature of the system, and hence expected to be case 
dependent. 
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